Magnetostructural effects and phase transition 
in Cr 2 03 under pressure 
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We have successfully calculated the electronic and structural properties of chromia (C^Oa) in 
the Local Spin Density Approximation (LSDA) . We predict a transformation from the corundum to 
the RI12O3 (II) structure around 15 GPa in the anti-ferromagnetic (AFM) phase as well as in the 
paramagnetic (PM) insulating state which occurs above the Neel temperature (TW). This transition 
is relevant to interpreting the optical anomalies observed in the absorption spectrum of ruby under 
pressure. We have modeled the structural properties of the PM state using a Landau-like expansion 
of the magnetostriction energy. This treatment correctly describes the structural anomalies across 
Tjv in the corundum phase and indicates that the AFM and PM insulating states should have 
distinct compressive behaviors. 
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Oxides of the 3d transition metals are a fascinating 
class of materials with amazingly diverse physical proper- 
ties. They have been a subject of intensive experimental 
and theoretical studies for many years jl]-[|]. A wide va- 
riety of computational techniques — density-functional 
theory with linear augmented plane wave method 
non-periodic embedded cluster approach periodic 
unrestricted Hartree-Fock method |]U — have been used 
to perform ab initio calculations of these compounds. 
In this paper we report a successful application of the 
first principles pseudo-potential plane- wave (PPPW) ap- 
proach |1(J to compute structural properties of and pre- 
dict a pressure induced structural transition in chromia 
(CV2O3), a typical anti-ferromagnetic insulator in this 
class of materials. There were a few experimental papers 
on chromia JTT|-p^| and only one ab initio calculation || . 
The periodic unrestricted LCAO Hartree - Fock method 
was used to calculate crystal parameters and elas- 
tic properties, as well as electronic structure of CV2O3. 
Their results are in a good agreement with experimen- 
tal data. The advantage of the PPPW approach is in 
the simplicity of the plane- wave basis set, which makes 
it easy to calculate ionic forces and lattice stresses. This 
allowed us to optimize dynamically the cell structure of 
chromia under pressure p"5| for various magnetic states. 

Our motivation to study CY2O3 is related to the high 
pressure behavior of ruby, i.e. Al 2 {i- x )Cr2 X 0^ (x < 
0.05). The pressure dependence of the fluorescence lines 
in ruby is widely used to determine pressure (the so-called 
ruby scale) in diamond-anvil-cell experiments. Alumina 
{AI2O3) and chromia exist in the corundum phase and 
form a completely isomorphous alloy system. Alumina 
has been shown to undergo a structural phase transition 
to the Rh 2 3 (II) phase around 80-100 GPa ||§hT§. A 
recent theoretical study JilJ of the effect of this structural 
transition on the optical spectrum of ruby indicated that 
the neighborhood of the chromium site, a color center, 



might be undergoing a severe distortion around 30 GPa. 
This hypothesis is suggested by the behavior of the op- 
tical absorption lines which display a small discontinuity 
at 30 GPa and resemble more closely the transitions pre- 
dicted in the high pressure Rh^Oj, (II) phase beyond 30 
GPa than those in the corundum phase. It was then 
anticipated that the cause of this distortion could be a 
similar phase transformation at lower pressures (< 30 
GPa) in chromia, the other end member of the alloy. 
At the moment there is no convincing experimental evi- 
dence for a phase transformation in CV2O3 in this pres- 
sure range pOfl . However, this may be due to the fact 
that the Rh^Oj, (II) phase has an X-ray diffraction pat- 
tern similar to corundum's Jl6| ], which makes it difficult 
to observe. 

In this paper we investigate this pressure induced 
transformation in chromia. However another interesting 
question arises: at room temperature chromia undergoes 
a change in magnetic phase under pressure. Its Neel tem- 
perature is T N = 308 K p) with a pressure dependence 
of dT N /dP = -16 K/GPa @. The paramagnetic (PM) 
state above T/v is also insulating; therefore, the effect of 
the magnetic transition on the structural properties are 
not expected to be dramatic. Nevertheless, structural 
anomalies around the Neel temperature are well known 
p3| , p4| and a realistic prediction of a possible phase tran- 
sition above 0.5 GPa should be carried out in the PM 
insulating phase with randomly oriented spins. Here we 
investigate from first principles this structural transition 
in the AFM, and in a hypothetical ferromagnetic (FM) 
phases. The structural properties of PM insulating state 
is then explored in relation to those of the AFM and 
FM phases using a phenomenological approach based on 
a Landau-like expansion of the magnetostriction energy. 
The predicted structural differences of this phase with re- 
spect to the AFM phase correlate well with the anomalies 
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Density of States (states/eV cell) 

FIG. 1. Calculated DOS for various magnetic phases of 
Cr 2 03 in the corundum structure. 

observed around the Neel transition and are very differ- 
ent from those of the PM metallic phase predicted by a 
standard LDA calculation. 

The crystal structure of CV2O3 at ambient conditions 
is corundum-like. It can be described as a hexagonal 
closed packed array of oxygens with two thirds of the 
octahedral sites filled with chromium atoms. The unit 
cell is rhombohedral and contains two formula units. In 
CY2O3 each chromium is left with three d-electrons, los- 
ing the other three to oxygens. The predominantly octa- 
hedral crystal field splits <i-orbitals into (approximately) 
a i2g-like triplet and an e 9 -like doublet. The lower triplet 
accommodates three electrons. The Cr 3+ local magnetic 
moment is 2.76/j,b, (2^] close to the spin only value of 
3/-J.B- The moments lie parallel to the c-axis, and are 
arranged in ferromagnetic layers perpendicular to the c- 
axis. The layers alternate up and down along the c-axis. 
In the insulating phase above Tjv chromia has randomly 
oriented spins. 

To investigate the effect of magnetism on the struc- 
tural properties we have performed three distinct calcu- 
lations: 1) a standard spin-polarized LSDA calculation 
p5f in the AFM phase of CriO-i in the corundum struc- 
ture; 2) same-type calculation in the FM phase, with the 
net magnetic moment of per Cr atom; 3) a standard 
non-spin-polarized LDA calculation in a PM phase. In all 
cases the lattice and internal degrees of freedom were dy- 
namically relaxed under pressure |l5| . The zero-pressure 
structures obtained correspond to various local minima 
of the LSD functional. 

The ground state at T = K is the AFM state with 
a band gap of » 1.5 eV (Fig. |l|-a) and a local magnetic 
moment of on chromium atoms (from straight band 
occupations). The overall band structure compares well 
with photo-emission data Namely, the C>2p and 

Ci^d band widths of 5 eV and 1 eV respectively, and C>2p- 
Cr 3( j band centers separation of 4 eV are in a good agree- 
ment with experimental values. However, as expected, 



TABLE I. Zero pressure structural parameters of chromia 
in the corundum structure (rhombohedral unit cell): lattice 
constant ao (A) and rhombohedral angle a (deg), internal 
atomic coordinates u{Cr) and u(O), bulk modulus Bo (GPa) 
and its pressure derivative B' . 



AFM LSDA FM LSDA PM LDA Experiment 



a 


5.366 


5.308 


5.688 


5.362 c 5.350 d 





55.17 


56.14 


47.23 


55.108 c 55.128 d 


u{Cr) 


0.347 


0.351 


0.337 


0.3475 c 0.3477 d 


u{0) 


0.557 


0.550 


0.583 


0.556 c 0.555 d 


B a 


251 ±6 


215 ±5 


300 ±8 


238 ±4 d 222 ±2 


B h 


261 ±2 


211 ± 5 


297 ±3 


231 ± 5 f 


B' b 


2.59 


2.73 


4.24 


2.0 ± l.l f 



a Second order finite strain equation of state (FSEoS) (B' =4). 

We used our data up to 15 GPa for this fitting. 

b Third order FSEoS was used with B' as free parameter. We 

include pressures up to 140 GPa to get a correct value of B' , 

while in Ref. Jl| the pressure range was not sufficient for a 

confident determination of B' . 

c Ref. H 

d Ref. m 

f Ref. fi|| 

the band gap is underestimated with respect to the ther- 
mal gap of 3.3 eV | 0T| . Zero-pressure equilibrium struc- 
tural parameters presented in Table | are also in good 
agreement with experimental data. The cohesive energy 
of E C oh = 6.1 eV/atom compares well with E co h — 5.55 
eV/atom from experiments. The FM state is found to be 
insulating as well with a band gap of w 0.9 eV (Fig.Q-b). 
After structural relaxation this state is only 35 meV/unit 
above the AFM ground state. Equilibrium lattice param- 
eters (Table |) are quite different from those in the AFM 
state, indicating a substantial influence of the magnetic 
state on the structural properties. The standard para- 
magnetic non- spin-polarized LDA calculation stabilizes 
CV2O3 in a metallic phase (Fig.Q-c), as expected, with 
structural properties considerably different from the ob- 
served ones (see Table |). This state is 2.25 eV/unit 
above the AFM ground state and cannot properly ac- 
count for the structural transition under consideration. 

Our description of the structural properties of the 
PM phase is based on the Landau-like expansion of the 
crystal deformation energy at a certain fixed pressure, 
E = XikjiUikUji + PikjiiiikMijMzi. Here u ik is a strain 
tensor, Mi, M2 are the magnetization vectors of the two 
AFM sublattices, and A and j3 are constant tensors. The 
first term represents pure elastic deformation, while the 
second describes the magnetostriction energy, i. e., the 
coupling of the magnetic and structural degrees of free- 
dom. Deformations conserving the corundum structure 
symmetry allow only for U — 1*33, uniaxial strain, and 
V = J2 u a> hydrostatic compression. We assume Mi, 
M2 remain parallel to the z-axis. The above expres- 
sion then simplifies to: E = X1V 2 + X 2 U 2 — X12VU + 
(PiV + P2U) ■ Mi z M 2z . The equilibrium configurations 
are such that U — uq ■ M\ Z M% Z , V — vq ■ M\ Z M2 Z , where 
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uq and v depend on A and j3. For the AFM and FM 
phases M\ z M 2z = — 1 and M\ z M 2z = 1 respectively. For 
the PM state the average product (M 1z M 2z )pm = 0. 
Therefore, the equilibrium lattice structural parameters 
of the PM phase are given by: Upm = {Uafm + Ufm)/2, 
Vp M = {Vafm + V FM )/2 where Uafm, U fm , Vafm, 
and Vfm have been determined from first principles. A 
similar procedure can be adopted for dealing with the 
RI12O3 (II) phase. In this case, the deformation energy 



is expressed as: E = ^ Xygu^us + ^2 Pi 



M lz M 2z 



where 7, 8 = xx, yy, zz. 

The predicted properties of the AFM and PM insu- 
lating states compare as follows: a) the difference in 
zero pressure lattice parameters (in the hexagonal cell 
description) between them, are similar to the anoma- 
lies observed around Tjy. Throughout the Neel transi- 
tion (AFM to PM), the calculated Aa H = +0.0131 and 
Acp = —0.11 A agree in sign and approximately in or- 
der of magnitude with experimental values of Aa# = 
+0.006A and Ac H = -0.0181 @. We believe this is 
evidence of the satisfactory description of the PM insu- 
lating state, b) The calculated compressive behavior of 
our PM insulating state compares well with the experi- 
mental behavior under pressure above T^ Jl2| , while our 
AFM calculation is in better agreement with the experi- 
ment of Lewis and Drickamer |ll| , which we suspect may 
have been carried out in the AFM phase. 

A summary of experimental data and our results is 
presented in Fig. ^. In Ref. [JllJ a substantial decrease of 
the rhombohedral angle with pressure was found (hexag- 
onal c-axis less compressive than a — b axes). In contrast, 
in Ref. jjj) a slight increase in this angle was observed 
(c-axis more compressive than a — b axes) . This discrep- 
ancy has been attributed to non-hydrostatic stresses in 
Ref. pT| . Our results suggest this discrepancy could be 
real if somehow in Ref. 1 1 chromia was kept in the AFM 
state. We predict here that the magnetic state affects no- 
ticeably the compressive behavior of chromia despite the 
insulating nature of both phases. Compression experi- 
ments well below and well above Tjy could help to clarify 
this situation. 

Now we deal with the structural transition. The high 
pressure Rh 2 0?, (II) phase has Pbna space group with an 
orthorhombic unit cell containing 20 atoms (four Cr203 
units). This phase is structurally similar to the corun- 
dum structure and may be described as having a different 
stacking of similar octahedral layers with a periodicity 
which is twice that of the corundum along the 'c' direc- 
tion |}2| . In the corundum structure the Cr06 octahedra 
share three edges while in the RI12O3 (II) structure they 
share only two. The relative stability of these structures 
in both magnetic states is shown in Fig. |^-a. We predict 
the corundum to Rh 2 0^ (II) transformation in chromia 
to take place at 14 GPa and 16 GPa in the AFM and 
in the PM insulating phases respectively, with fractional 
volume changes at the transition of rj —2%. This tran- 
sition would also take place in the FM phase if it were 
somehow stabilized. 
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FIG. 2. Pressure dependence of the CY2O3 rhombohedral 
cell constant a and angle a, calculated for AFM, FM and 
PM phases compared to experimental measurements from 
Ref. jO} (Lewis 66) and Ref. @ (Sato79). 



Fig. g-b displays the average radius of the first co- 
ordination shell around chromium in AFM chromia and 
in ruby (from Ref. jL9)). The average Cr-0 bond-length 
in chromia increases across the corundum to RI12O3 (II) 
transition by an amount similar to that required to ex- 
plain the optical anomalies in ruby under pressure, which 
can be explained by a decrease in crystal field splitting. 
This verifies that this presumable rearrangement could 
arise from a preference of chromia for the RI12O3 (II) 
phase above m 15 GPa. The structural constraint im- 
posed by the alumina host structure should naturally 
hinder the atomic rearrangement around the color cen- 
ters until higher pressures, for instance 30 GPa. 

These results should stimulate further experimental 
and theoretical work. The prediction of distinct com- 
pressive behaviors in chromia above and below Tjy and 
the structural phase transformation near 15 GPa await 
experimental confirmations. The latter, if verified, makes 
ruby an interesting study case: an isomorphous alloy in 
which both end members undergo the same structural 
transition but at very different pressures. Intermediate 
compositions should undergo similar transitions at inter- 
mediate pressures. However, before the transformation 
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(a) Pressure dependence of the enthalpy for 
RI12O3 (II) structure, relative to the corundum in AFM and 
PM phases; (b) Average radii of the first coordination shell 
around chromium in chromia and ruby Jl!| , and aluminum in 
alumina 



manifests macroscopically it could be nucleating around 
one of the components, even in the impurity limit. The 
possibility of investigating this phenomenon in ruby by 
EXAFS or anomalous X-ray scattering is fascinating. 
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